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摘 要: 针对 当前 麦克 风 位 置 估计 方法 出 现 较 少 且 只 对 平面 的 麦克 风 阵 列 进行 位 置 估计 的 问题 ， 提 出 了 一 种 基于 上 声 
源 阵列 的 空间 麦克 风 定 位 方法 。 该 方法 将 四 个 声 源 按照 正四 面体 形式 排列 , 在 利用 其 对 称 性 保证 定位 全 向 性 的 同时 ， 
根据 时 延 估计 的 测 距 方法 计算 麦克 风 与 各 个 声 源 之 间 的 距离 ， 之 后 根据 几何 关系 确定 麦克 风 位 置 坐标 ， 实 现在 线 的 
基于 声 源 阵列 的 空间 麦克 风 位 置 估计 。 通 过 进行 MATLAB 仿真 实验 ， 验 证 了 该 方法 的 可 行 性 和 有 效 性 
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Microphone spatial positioning method based on sound source array 


Chen Xiaohui, Sun Hao', Zhang Heng, Zhai Baoshuo 
(School of Artificial Intelligence, Hebei University of Technology, Tianjin 300130, China) 


Abstract: Currently, there are fewer microphone position estimation methods and only position estimates are made for a 
planar microphone array. Aiming at this problem, this paper proposed a spatial microphone positioning method based on 
sound source array. The method arranges four sound sources in the form of a regular tetrahedron. While using its symmetry 
to ensure the omnidirectionality of the positioning, the distance between the microphone and each sound source is calculated 
according to the ranging method of the time delay estimation. The microphone position coordinates are then determined 
based on the geometric relationship to realize the spatial microphone position estimation based on the sound source array for 
online positioning. In this paper, MATLAB simulation experiment is carried out, and the experimental results verifies the 
feasibility and effectiveness of the proposed method. 
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融合 得 到 麦克 风 位 置 。 但 该 方法 的 信号 采集 操作 较 繁 琐 ， 效 


0 引言 率 较 低 ， 只 外 离线 计算 ， 无 法 实现 麦克 风 位 置 的 在 线 估计 。 
伴随 着 科学 技术 日 新 月 异 的 发 展 ， 传 感 器 技术 的 应 用 领 “文献 [12] 通 过 选择 三 个 不 同 频带 的 声音 信号 作为 三 个 声 源 的 
域 逐步 扩大 ， 其 所 面 对 的 工作 环境 也 愈 来 愈 复杂 巾 。 对 麦克 ”发 声 信 号 ， 在 麦克 风 位 置 估计 过 程 中 一 次 性 采集 并 分 离 出 三 
7 置 进行 准确 的 定位 能 够 帮助 分 布 式 麦 克 风 阵列 信号 处 理 ， 个 声 源 ， 实 现 麦 克 风 位 置 的 在 线 估 计 。 将 基于 能 量 和 时 延 的 

过 克 风 阵列 在 语音 增强 、 声 源 定位 、 声 源 分 离 等 应 用 中 ” 测 距 方 法 结合 使 用 ， 先 后 进行 距离 的 粗 、 精 估计 ; 采用 三 


起 着 重要 作用 。 针 对 如 何 对 麦克 风 阵 列 信 号 个 7 进行 处 理 的 间 ” 形 质 心算 法 获取 阵列 中 各 麦克 风 的 位 置 ， 对 各 种 无 解 情况 均 

究 虽 已 近 十 年 ， 但 大 多 集中 于 国外 ， 并 且 由 于 刚 起 步 不 有 相应 的 处 理 ， 该 方法 保证 了 测 距 具 有 较 小 的 计算 量 和 较 高 
还 仍 处 于 探索 阶段 。 然 而 在 国内 ， 对 于 该 问题 的 研究 则 ”的 精度 ， 同 时 保证 定位 结果 的 持续 性 和 稳定 性 。 但 此 方法 只 
少 。 因 此 最 终 研究 成 果 不 多 。 同 时 ， 准 确 得 到 麦克 风 位 。 针对 平面 的 麦克 风 阵 列 进 行 了 位 置 估计 ， 无 法 实现 空间 分 布 
E 标 可 以 帮助 解决 当前 机 器 人 领域 众多 问题 ， 例 如 ， 当 机 ”麦克 风 的 位 置 估 计 。 
人 处 于 未 知 的 环境 中 ， 对 自身 位 姿 未 知 时 ， 如 何 进行 自 定 为 解决 现 有 技术 存在 的 上 述 问 题 ， 本 文 提 出 了 一 种 精度 
位 和 环境 识别 ， 即 同时 自 定位 和 地 图 构建 由 一 直 都 是 智能 机 较 高 、 计 算 复 杂 度 较 小 、 并 能 够 实现 在 线 定位 的 基于 声 源 阵 
器 人 技术 研究 的 热点 和 难点 外。 然而 从 机 器 人 听觉 系统 的 角 ” 列 的 空间 麦克 风 定 位 方法 。 在 一 定 的 空间 范围 内 ， 随 机 分 布 
度 出 发 ， 确 定 出 待定 位 麦克 风 的 位 置 坐标 ， 将 机 器 人 技术 与 ” 一 定数 量 的 麦克 风 ; 声 源 分 布 单元 为 四 个 声 源 按照 四 元 三 维 
机 器 人 所 处 环境 中 的 分 布 的 麦克 风 位 置 进行 有 机 结合 ， 可 帮 。 ”阵列 ， 特 别 是 正四 面体 的 方式 排列 分 布 ， 声 源 分 布 单元 依次 
助 解决 机 器 人 同时 自 定位 与 构建 环境 地 图 的 问题 ， 从 而 实现 。 产生 声音 信号 ， 接 收 到 四 个 声音 信号 的 麦克 风 定 义 为 待定 位 
机 器 人 全 自主 移动 等 "4。 因 此 ， 确 定 麦克 风 阵 列 中 的 麦 殉 风 ”麦克 风 ， 之 后 采用 时 延 估 计 法 获取 参考 麦克 风 与 待定 位 麦克 
位 置 是 一 个 必须 解决 的 技术 难题 。 风 接 收 到 声音 信号 的 时 间 差 ( 声 达 时 间 差 ) 信息 ， 最 后 根据 
文献 [11] 中 采用 三 个 不 同位 置 的 校正 声 源 来 校正 麦克 风 几何 关系 确定 各 个 待定 位 麦克 风 位 置 。 此 方法 利用 正四 面体 
阵列 中 的 麦克 风 位 置 。 其 方法 是 将 一 个 声 源 先后 放置 于 三 个 ” 的 对 称 性 ， 保 证 了 定位 全 向 性 ， 同 时 保证 信号 被 麦克 风 均 匀 
位 置 ， 分 三 次 进行 信号 的 采集 与 计算 ， 最 后 将 三 次 计算 结果 ”接收 ,从 而 满足 定位 精度 ,保证 定位 精度 的 持续 性 和 稳定 性 。 
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维 空间 中 的 麦克 风 的 在 线 定位 。 


了 
1 ” 声 源 阵列 模型 
个 


布 


本 文 将 声 源 分 布 单元 按照 四 元 三 维 阵列 分 布 ， 即 构建 一 


掉 体 的 形式 空间 分 


三 维 空间 坐标 系 ， 将 四 个 声 源 按照 正 


。 可 利用 其 对 称 性 保证 信号 被 麦克 


3 


均匀 接收 ， 从 而 满足 


4E 


为 


基准 。 声 源 阵列 模型 如 图 1 所 示 。 


立 精 度 ， 同 时 空间 分 布 的 麦克 风 位 置 表示 都 将 以 此 坐标 系 


图 1 声 源 阵列 模型 


Fig.1 Themodel of sound source array 


图 中 ， 在 以 o 为 原点 的 三 维 空间 坐标 系 (x、y、z) 中 ， 


S1、S2、S3、S4 分 别 为 四 个 声 源 ，C0 为 参考 麦克 风 ， 其 所 
处 位 置 距离 各 个 声 源 小 于 10 米 ，M1 设 定 为 第 一 个 待定 位 麦 


克 风 ， 即 所 要 进行 定位 的 目标 待定 位 麦克 风 之 一 ,Ll 为 参考 


麦克 风 C0 到 第 一 个 声 源 S1 的 距离 ,已 知 且 固定 不 变 , 同 理 ， 


参考 麦克 风 C0 到 其 他 声 源 (S2、S3、S4) 的 距离 分 别 为 L2、 
L3、L4 也 已 知 且 固定 不 变 ，d11 为 第 一 个 待定 位 麦克 风 M1 
到 第 一 个 声 源 S1 的 距离 ， 它 是 由 S1 到 达 C0 与 M1 的 时 间 


差 乘 以 声速 计算 得 到 的 距离 差 ， 然 后 再 根据 已 知 距离 L1 所 
求 得 的 ， 依 此 类 推 ， 可 计算 出 M1 到 其 他 声 源 的 距离 ， 这 是 


定 


2 


器 所 接收 信号 的 时 间 差 进行 估计 ， 来 确定 
新 源 的 距离 、 方 位 、 速度 和 移动 方向 等 3 。 
[将 时 延 估计 方法 分 为 如 下 几 类 :相位 法 、 


条 件 和 信号 的 特性 可 


出 M1 空间 位 置 的 必要 条 件 。 
时 延 估 计算 法 


时 延 估 计 主 要 指 利用 信号 处 理 的 理论 和 方法 对 


了 
其 他 相关 参量 ， 如 


按照 测量 环境 、 


双 谱 法 、 相 关 法 、 自 适应 滤波 器 参数 模型 法 等 。 甚 中， 相关 


声 源 分 布 单元 依次 产生 声音 信号 ， 参 考 麦 克 风 与 待定 位 
麦克 风 接 收 到 声音 信号 。 这 里 假设 麦克 风 接 收 到 的 声音 信号 
只 折 


模型 是 理想 模型 ， 即 上 只 考虑 环境 噪声 ， 这 里 环境 噪 用 高 其 
噪声 近似 代替 ， 且 假设 噪声 与 噪声 之 间 、 噪 声 与 声 源 之 间 


白 


法 是 最 经 典 的 且 常 用 的 一 种 时 延 估 计 方 法 。 
2.1 基本 互相 关 法 
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本 文 利用 声 达 时 间 差 数据 ， 通 过 MATLAB 仿真 实验 ， 实 现 如 式 〈2) 所 示 。 


Rs (7)=E[x (x(t—7)] (2) 
将 式 (1) 代入 式 (2) 可 得 


R,, (7) 过 aR, (z Tt ) 二 CR (7) 
+ QR,, (tT—Ty)+ R,, (7) 


G3) 


其 中 : RR,(7) 是 s(t) 的 自 相关 函数 ，Rw 是 m0 与 (0 的 互相 


关 函 数 ，&w 和 Rs 是 噪声 与 声 源 的 互相 关 函 数 ，w 是 第 i、 j 


个 麦克 风 接 收 到 声音 信号 的 相对 时 间 延 迟 。 因 为 噪声 与 噪声 
过 间 、 噪 声 与 声 源 之 间 互 不 相关 ， 则 式 〈3) 可 化 简 为 


R,, (7)= 00)R, (tr-zt;) (4) 


根据 自 相关 函数 的 性 质 可 知 当 T-5=0 时 ，RsstW) 达 到 最 


大 值 0。 则 时 延 估 计 不 便 可 由 两 路 麦克 风 接 收 信号 的 互相 关 
函数 和 其 最 大 值 位 置 所 在 点 计算 而 得 。 
2.2 广义 互相 关 法 
广义 互相 关 算 法 是 目前 应 用 最 广泛 的 一 种 时 间 估 计 方 法 
09。 针 对 实际 环境 中 噪声 与 混 响 大 量 存 在 ， 进 而 造成 互相 关 
函数 最 大 峰值 不 明显 或 者 峰值 出 现 多 个 的 情况 ， 广 义 互相 关 
以 其 抑制 噪声 与 混 响 的 能 力 相 对 基本 互相 关 算 法 在 峰值 检测 
难 易 程度 方面 性 能 较 好 。 对 式 (4) 进行 傅 里 叶 变 换 可 得 


pis (@) = 000s (oO)e ™™ + 9 (OO) (5) 


其 中 : Vs Co) ，96w (8) ，Vm (89) 分 别 为 R(T)，Rss (7) 与 Ri (TD) 


所 对 应 的 功率 谱 , 其 中 ，&w (为 n(n) 和 (7) 的 互相 关 函 数 。 


进而 ,， 对 式 (5) 进行 加 权 后 做 道 变换 ， 则 可 得 到 广义 互相 关 
函数 如 式 (6) 所 示 。 


已 (= jos (0)0, (0)emdo (6) 


其 中 : 8(%@) 为 广义 加 权 函 数 17。 加 权 函 数 有 Roth、SCOT 
(平滑 相干 变换 )、PATH (相位 变换 ) 等 。 由 于 它们 有 着 不 
同 的 特性 ， 应 根据 不 同 应 用 场景 选择 不 同 的 加 权 函 数 。 由 于 
广义 互相 关 以 其 抑制 噪声 与 混 响 的 能 力 相 对 基本 互相 关 算 法 
在 峰值 检测 难 易 程度 方面 性 能 较 好 ， 因 此 本 文选 取 广 义 互相 
关 时 延 估 计算 法 。 并 且 PATH 加 权 函 数 在 信 噪 比 不 断 变化 的 
情况 下 有 着 较 小 的 波动 性 ， 抗 噪声 性 能 也 较 好 08， 因 此 其 中 
应 用 PATH 加 权 函 数 进行 定位 仿真 分 析 。 
3 麦克风 位 置 估计 

在 实际 环境 中 ， 将 一 定数 量 的 麦克 风 准 确 得 摆 放 在 同一 
平面 内 有 一 定 的 难度 ， 因 此 对 同一 平面 内 的 麦克 风 位 置 估计 
将 会 产生 一 定 的 误差 。 而 将 麦克 风 在 同一 空间 内 随机 分 布 相 
对 容易 ， 即 降低 了 设备 组 装 等 的 困难 程度 ， 并 且 实 现 麦克 风 


闻 位 置 估计 可 帮助 解决 传感器 技术 、 机 器 人 技术 等 领域 众 


互 不 相关 。 麦 克 风 接收 到 信号 的 数学 模型 如 式 (1)〉 所 示 。 


X(t)=0s(t—T)+n (1) 


DH. 


™ 


ti、 


数 、 


(1) 


: 声 源 信号 为 s(n) ,第 i 个 麦克 风 接 收 到 的 信号 为 *(7) ,%、 


ni(?) 分 别 表示 声 源 信号 到 达 第 i 个 麦克 风 的 幅度 衰减 系 
时 间 延 迟 和 环境 噪声 。 第 i 和 j 个 麦克 风 的 互相 关 函 数 


多 问题 。 
本 文 所 述 的 麦克 风 空 间 位 置 估计 系统 包括 声 源 分 布 单元 、 
考 麦 克 风 、 待 定位 麦克 风 以 及 麦克 风 定 位 单元 。 在 一 定 的 
范围 内 ， 随 机 分 布 一 定数 量 的 麦克 风 ; 声 源 分 布 单元 为 
声 源 按照 第 一 节 所 述 形式 排列 ， 并 依次 发 声 ， 共 产生 四 
音信 号 ; 将 在 空间 中 分 布 的 所 有 麦克 风 中 能 接收 到 四 个 


他 马凯 迷 
相依 丁 
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声音 信号 的 麦克 风 定 义 为 待定 位 麦克 风 ， 即 该 四 个 声音 信号 
被 待定 位 麦克 风 和 参考 麦克 风 所 接收 ， 参考 麦克 风 与 待定 位 
麦克 风 接 收 声 音信 号 的 时 间 差 ( 声 达 时 间 差 ， 传 递 给 麦克 风 
定位 单元 ， 经 过 上 节 所 述 广 义 互 相关 时 延 估 计 法 获取 声 达 时 
间 差 之 后 ， 麦 克 风 定位 单元 根据 声 达 时 间 差 以 及 参考 麦克 风 
与 待定 位 麦克 风 之 间 的 位 置 关 系 确定 出 待定 位 麦克 风 的 位 置 。 
基于 声 源 阵列 的 空间 麦克 风 定 位 方法 流程 图 如 图 2 所 示 。 


S2 ls 人 > TDOAL 
Oo 时 延 
i 2 >» TDOA2 二 网 定位 
@ ee | 败 估 计 > i 
a 法 | FLmon 站 元 
54 
® >| 入 雪 充 | 
RN >» TDOA4 


图 2 基于 声 源 阵列 的 空间 麦克 风 定 位 方法 流程 图 
Fig.2 Flow chart of spatial microphone positioning method based on 


sound source array 


如 图 2 所 示 ，S1、S2、S3、S4 分 别 为 产生 声音 信号 的 四 图 3 ”四 个 球 相交 于 一 点 示意 图 
个 声 源 ， 产 生 的 声音 信号 被 参考 麦克 风 和 待定 位 麦克 风 1、 Fig.3 Schematic diagram of four balls intersecting at one point 
待定 位 麦克 风 2… 待 定位 麦克 风 V 所 接收 ， 利 用 经 典 的 时 延 图 4 为 四 个 球 相 交 于 一 点 在 三 维 空间 中 较为 直观 的 表示 ， 


AAA a 

估计 法 中 的 广义 互相 关 函 数 法 来 估计 四 个 声 源 依次 产生 的 j 以 其 中 一 个 待定 位 麦克 风 为 例 。 图 中 红色 圆 点 即 为 某 一 个 待 
音信 号 s(n)(i=12,3,4) 到 达 参 考 麦 克 风 与 待定 位 麦克 风 的 时 ”定位 麦克 风 的 位 置 。 
间 差 即 声 达 时 间 差 , 即 为 图 中 的 TDOA1、TDOA2、TDOA3、 
TDOA4, 其 中 ,上 表示 各 个 声 源 发 声 对 应 的 时 间 ， 到 12,3,4 表 
示 四 个 声 源 ， 进 行 “麦克 风 定 位 单元 ”的 操作 ,将 TDOAI1、 
TDOA2、TDOA3、TDOA4 作为 输入 信号 ， 输 出 待定 位 麦克 
风 在 三 维 空间 中 的 坐标 ， 包 括 待定 位 麦克 风 1 坐标 、 待 定位 
麦克 风 2 坐标 … 待 定位 麦克 风 N 坐标 ， 由 此 完成 基于 声 源 阵 
列 的 空间 麦克 风 定 位 ,麦克 风 定 位 单元 的 工作 流程 描述 如 下 : 

定义 AT(i=1,2..….N;j=1,2,3,49 为 参考 麦克 风 与 第 i 个 待定 
位 麦克 风 接 收 到 第 j 个 声音 信号 的 时 间 差 ， 定 义 声音 在 空气 
中 传播 速度 为 C， 将 信号 到 达 时 间 差 乘 以 声速 C， 得 到 第 j 
个 声 源 与 参考 麦克 风 的 距离 和 与 第 i 个 待定 位 麦克 风 的 距离 
的 差 值 如 式 〈7) 所 示 。 


AD =AT; xC (7) 


根据 参考 麦克 风 到 第 j 个 声 源 的 距离 50=142,3,4， 以 及 
第 j 个 声 源 到 参考 麦克 风 与 第 i 个 待定 位 麦克 风 的 距离 差 Am ， 
获得 第 i 个 待定 位 麦克 风 到 第 j 个 声 源 的 距离 必 

djy=L;+AL;=L,+AT;xC (8) 

根据 几何 关系 ， 通 过 求解 分 别 以 四 个 声 源 位 置 为 球 心 ， 
以 待定 位 麦克 风 到 相对 应 的 声 源 的 距离 为 半径 的 四 个 球 的 交 图 4 四 个 球 相交 于 一 点 立体 示意 图 
点 ， 来 确定 第 i 个 麦克 风 在 空间 环境 中 的 位 置 。 即 Fig.4 Stereoscopic diagram of four balls intersecting at one point 

(x—Xn) +(y— +(z—a) = dr 4 ”实验 与 分 析 


(Xx—X) +(y—y) +(z— 2) = 4d (9) 

xn) +(y-y) +(z—z) = ds 为 了 验证 方法 的 可 行 性 和 有 效 性 ， 基 于 上 述 广义 互相 关 

C0 FY) PE dy 时 延 估计 法 获取 的 声 达 时 间 差 数据 ， 本 文 在 MATLAB 仿真 

于 四 个 声 源 位 置 在 空间 中 以 正四 面体 形式 分 布 ， 因 而 环境 下 进行 了 仿真 实验 ,进行 了 如 下 仿真 测试 。 设 定 一 个 长 、 

由 声 源 两 两 构成 的 向 量 线性 无 关 ， 从 而 保证 上 述 方程 只 有 一 宽 、 高 分 别 为 60 m、60 m、20 m 的 空间 区 域 ， 如 图 5 中 人 举 
解 


配 


个 解 ， 该 解 即 为 四 个 球 的 交点 ， 此 交点 坐标 即 为 第 i 个 待定 “， 标 轴 所 示 ， 在 该 空间 区 域内 ， 图 中 四 个 空心 圆 口 示 为 四 个 声 
位 麦克 风 在 空间 中 的 位 置 坐标 .图 3 为 其 几何 模型 解释 :S1、 源 的 坐标 ， 位 置 分 别 为 (6,0,0) m、 (-3.3V3.0) m、 (-3,-3V3,0) m 
S2、S3、S4 分 别 为 四 个 声 源 ， 以 第 一 个 待定 位 麦克 风 M1 为 与 (0,0,6V2)m; 参考 麦克 风 友 坐标 为 (0,0,0) m， 即 选择 坐标 原 
例 ，dl1、d12、d13、d14 分 别 为 M1 到 四 个 声 源 的 距离 。 以 ” 点 为 参考 麦克 风 友 的 位 置 ， 图 中 空间 区 域 中 有 200 个 随机 分 
声 源 位 置 为 球 心 ， 相 对 应 的 待定 位 麦克 风 到 声 源 的 距离 为 半  ” 布 的 麦克 风 十 ， 经 验证 ， 其 中 有 两 个 接收 到 四 个 声音 信号 的 
径 作 球 ， 由 于 声 源 在 空间 中 以 正四 面体 形式 分 布 ， 因 而 由 j 麦克 风 ， 即 有 两 个 待定 位 麦克 风 @， 因 其 是 200 个 随机 分 布 

的 

的 


[之 


两 两 构成 的 向 量 线性 无 关 ， 从 而 该 四 个 球 在 空间 中 只 有 一 麦克 风 十 中 的 两 个 ， 所 以 在 图 中 显示 两 个 待定 位 麦克 凡 @ 
交点 ， 此 交点 即 为 待定 位 麦克 风 M1 的 坐标 ， 如 图 黑色 三 图 形 “@ ”肯定 是 与 200 个 代表 随机 分 布 的 麦克 风 十 中 的 
攻 所 示 。 两 个 “十 ”图 形 重 


下 全 东 


| 内 


上 孝 


， 重 肢 后 的 图 形 在 该 图 5 中 有 显示 。 根 
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据 本 文 方法 对 这 200 个 随机 分 布 的 麦克 风 中 的 两 个 待定 位 麦 
克 风 的 位 置 进行 了 估计 ， 估 计 的 结果 如 该 图 5 中 的 “ 口 ” 所 
示 ,“ 口 ?所 在 的 位 置 就 是 本 文 方法 对 麦克 风 估 计 的 位 置 ， 换 
句 话说 ,“@” 所 在 的 位 置 是 需要 对 其 进行 定位 的 也 就 是 待定 
位 麦克 风 的 实际 位 置 。 由 基于 声 源 阵列 的 空间 麦克 风 定 位 方 
个 待定 位 麦克 风 位 置 估计 出 来 ， 声 源 位 置 、 参 考 麦 
克 风 、 随 机 分 布 的 麦克 风 、 待 定位 麦克 风 以 及 待定 位 麦克 风 
位 置 估 计 的 分 布 如 图 5 所 示 。 


十 随机 分 布 的 雪 克 风 
交 参考 支 克 风 
a 日 待定 位 麦克 风 
宝生 待定 位 麦克 风 位 置 估计 
10 + + 〇 ”下 源 他 
人 b + 
+ + + 让 中 
+ 十 四 ++ + + 
We 二 HH 让 村 检 + 十 
~ 出 下 十 + 二 有 由 
类 + 村 :+ + 十 再 
2 0 + + + Or + 十 十 本 
吧 上 二 + 再 + 二 加 
十 十 十 十 
5 + 十 4 
证 
+ 


宽 ( 米 ) 30. .0 长 ( 米 ) 


图 5 仿真 结果 
Fig.5 Simulation result 
空间 麦克 风 定 位 结果 如 表 1 所 示 。 
表 1 空间 麦克 风 定 位 结果 分 析 


Table 1 Analysis of spatial microphone positioning results 
序号 克 风 坐标 真 值 /ma ” 信 品 比 /dB 麦克 风 坐 标 测量 值 /m 误差 /m 
1 (2.316,-0.744,5.675 ) 10 (2.307,-0.741,5.666) 0.012 
2 (1.313,-2.844,-0.205 ) 10 (1.320,-2.851,-0.221) 0.010 


从 表 1 可 以 看 出 ， 考 虑 在 一 定 的 信号 信 品 比 下 ， 麦 克 风 
在 三 维 空间 坐标 系 下 的 位 置 坐标 都 能 够 很 好 地 被 估计 出 来 ， 


定位 误差 较 小 。 


5 ”结束 语 

本 文 针 对 麦克 风 位 置 估计 问题 ， 提 出 了 一 种 基于 声 源 阵 
列 的 空间 麦克 风 定 位 方法 。 该 方法 采取 四 个 声 源 按照 正四 面 
本 形式 排列 的 方式 作为 声 源 分 布 单 元 ， 根 据 广义 互相 关 时 延 
占 计 方 法 获得 声音 信号 到 达 参 考 麦 克 风 与 待定 位 麦克 风 的 时 
同 差 ， 利 用 声 达 时 间 差 定 的 运算 关系 的 得 出 待定 位 麦克 
风 与 各 个 声 源 之 间 的 距离 ， 之 后 根据 几何 关系 ， 通 过 求解 分 
别 以 四 个 声 源 位 置 为 球 心 ， 以 待定 位 麦克 风 到 相对 应 的 声 源 
的 距离 为 半径 的 四 个 球 的 交点 ， 来 确定 该 麦克 风 在 空间 环境 
中 的 位 置 。 由 于 四 个 声 源 为 对 称 性 分 布 ， 因 而 保证 了 定位 全 
向 性 , 同时 保证 信号 被 麦克 风 均 匀 接 收 ,从 而 满足 定位 精度 ， 
保证 定位 精度 的 持续 性 和 稳定 性 。 此 方法 精度 较 高 、 计 算 复 
杂 度 较 小 、 并 能 够 实现 在 线 定 位 的 基于 声 源 阵列 的 空间 麦克 
风 位 置 估 计 。 对 麦克 风 进 行 空间 层面 的 定位 是 一 种 新 颖 且 极 
潜力 的 新 技术 ， 如 何 将 该 技术 应 用 到 声音 以 及 机 器 人 等 领 
域 是 一 个 值得 深入 研究 的 问题 。 


参考 文献 : 


[1] 和 孙 昊 . 基于 声 达 时 间 差 的 移动 机 器 人 上 声 源 目 标定 位 方法 研究 [D]. 
天 津 : 河北 工业 大 学 , 2012. (Sun Hao. Study on sound source 


localization for mobile robot based on time delay of arrival [D]. 
Tianjin:Hebei University of Technology, 2012. ) 
[2] 陈 国 良 , 徐 扬 , 黄 晓 琴 . 基于 正四 棱锥 结构 的 机 器 人 声 源 定位 系统 


六 Rinyay nA 化 革 日 于 || 
ChinaXiv 合 作 期 刊 


第 37 卷 第 5 期 


研究 [加 . 计算 机 应 用 研究 , 2019,36(5) . (Chen Guoliang, Xu Yang, 
Huang Xiaoqin. Sound source localization System for robot based on 
rectangular pyramid structure [J]. Application Research of Computer, 
2019,36 (5).) 

[3] Yook D, Lee T, Cho Y. Fast sound source localization using two-level 
search space clustering [J]. IEEE Trans on Cybern, 2017, 46(1): 20-26. 

[4] Li X, Liu H. Sound source localization for hri using FOC-based time 
difference feature and spatial grid matching [J].IEEE Trans on 
Cybernetics, 2013, 43(4): 1199-1212. 

[5] Benesty J, Chen J, Huang Y, et al. On microphone-array beamforming 
位 om a MIMO acoustic signal processing perspective [J]. IEEE Trans on 
Audio Speech & Language Processing, 2007, 15(3): 1053-1065. 

[6] Silva B D, An B, Steenhaut K, et al. Design considerations when 
accelerating an FPGA-based digital microphone array for sound-source 
localization [J]. Journal of Sensors, 2017 (4): 1-20. 

[7] 张 雷 项 , 张 兴 敢 , 刘 超 . 麦克 风 阵 列 声 源 定位 中 时 延 估计 的 改进 [J. 
南京 大 学 学 报 : 自然 科学 版 ，2015, 51(1): 25-30. (Zhang Leiyue， 
Zhang Xinggan, Liu Chao. The improvement of time delay estimation 
in the microphone array Sound localization System [可 .Journal of 
Nanjing University: Natural Sciences, 2015, 51(1): 25-30. ) 

[8] 柯 文 德 , 蔡 则 苏 , 李 家 兰 . 机 器 人 同时 定位 与 地 图 构建 技术 研究 [中 . 
计算 机 应 用 研究 , 2010, 27(4): 1216-1219. (Ke Wende, Cai Zesu, Li 
Jialan. Research of simulation localization and mapping in robot [J]. 
Application Research of Computers, 2010, 27 (4): 1216-1219. ) 

[9] 罗 元 ， 苏 琴 ， 张 发 ,等 . 基于 优化 RBPF 的 同时 定位 与 地 图 构建 [J]. 
华中 科技 大 学 学 报 : 自 然 科 学 版 , 2016, 44(5): 30-34. (Luo Yuan, Su 
Qin, Zhang Yi, et al. Simultaneous localization and mapping 
implementation based on optimized RBPF [J]. Journal of Huazhong 
University of Science and Technology: Natural Science Edition, 2016, 
44(5): 30-34. ) 

[10] Evers C, Naylor P A. Acoustic SLAM [J]. IEEE/ACM Trans on Audio 
Speech & Language Processing, 2018, PP(99): 1-1. 

[11] 肖 华 . 麦克 风 阵列 的 校正 方法 研究 [D]. 成 都 : 电子 科技 大 学 , 2008. 
(Xiao Hua. Research on correction method of microphone array [D]. 
Chengdu:University of Electronic Science and Technology of China, 
2008. ) 

[12] 王 舒 文 . 分 布 式 麦克 风 阵 列 定位 方法 研究 [D]. 大 连 :大 连理 工大 学 ， 
2013.(Wang Shuwen. Study on distributed microphone array 
localization method [D]. Dalian:Dalian Uuniversity of Technology, 
2013) 

[13] 行 鸿 产 ， 唐 娟 ,时 延 估 计 方 法 的 分 析 [四 . 声学 技术 ，2008，27(]): 
110-114. (Hang Hongyan, Tang Juan. Analysis and Survey of algorithms 
for time-delay estimation [J]. Technical Acoustics, 2008, 27(1): 
110-114. ) 

[14] Chen Z, Li Z, Wang S, et al. A microphone position calibration method 
based on combination of acoustic energy decay model and TDOA for 
distributed microphone array [J]. Applied Acoustics, 2015, 95: 13-19. 

[15] 行 鸿 关 ， 赵 守 国 ， 拨 继 征 , 等 . 广义 相关 时 延 估计 算法 的 自 适 应 实现 
形式 [J]. 西安 石油 大 学 学 报 : 自 然 科 学 版 , 2001, 16 (6): 47-50. (Hang 
Hongyan, Zhao Shouguo, Di Jizheng, et al. Adaptive realizing form of 
generalized correlation algorithm for time delay estimation [J]. Journal 
of Xi’” an Petroleum Institute: Natural Science Edition, 2001, 16(6): 
47-50. ) 

[16] 茅 惠 达 ， 张 玲 华 . 声 源 定 位 中 广义 互相 关 时 延 估 计算 法 的 研究 [J]. 
计算 机 工程 与 应 用 ，2016，52(22): 138-142. (Mao Huida，Zhang 


Linghua. Research on generalized cross correlation algorithm for time 


201904.00036v1 


chinaXiv 


delay estimation in sound Source localization [J]. Computer 
Engineering and Applications, 2016, 52(22): 138-142. ) 

[17] 张 莉 . 基于 传声器 阵列 的 声 源 定位 方法 研究 [D]. 成 都 : 电子 科技 大 
学 , 2007. (Zhang Li. Research on sound source localization method 
based on microphone array [D]. Chengdu:University of Electronic 


Science and Technology of China, 2007. ) 


ChinaXiv 合 作 期 刊 


录用 定稿 陈晓辉 ， 等 : 基于 声 源 阵列 的 空间 麦克 风 定 位 方法 研究 第 37 卷 第 5 期 


[18] 孙 洋 ， 徐 慨 ， 杨 海 亮 . 基于 广义 互相 关 时 延 估 计算 法 的 性 能 分 析 [. 
计算 机 与 数字 工程 ，2013, 41(1): 33-34. (Sun Yang，Xu Kai，Yang 
Hailiang. Performance analysis of time-delay estimation based on 
generalized cross-correlation algorithm [J]. Computer and Digital 
Engineering, 2013, 41(1): 33-34. ) 


